The effect of an external magnetic field on the gas-liquid transition in the Ising spin 

fluid 
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The theoretical phase diagrams of the magnetic (Ising) lattice fluid in an external magnetic 
field is presented. It is shown that, depending on the strength of the nonmagnetic interaction 
between particles, various effects of external field on the Ising fluid take place. In particular, at 
moderate values of the nonmagnetic attraction the field effect on the gas-liquid critical temperature 
is nonmonotoneous. A justification of such behavior is given. If short-range correlations are taken 
into account (within a cluster approach), the Curie temperature also depend on the nonmagnetic 
interaction. 
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Anisotropic liquids are very sensitive matters. Such 
are nematic liquid crystals and ferrofiuids. Many efforts 
have been made in order to investigate effects of shape 
and flexibility of the molecules, of long and short-range 
interactions on the properties of anisotropic liquids. Ex- 
ternal field effects are still worthier of attention, because 
the application of an external field allows to change prop- 
erties of anisotropic liquids dynamically (in contrast to 
the effects of molecules' shape, etc., which are static). 

In ferrofiuids an external magnetic field removes the 
magnetic order-disorder transition, nevertheless the first 
order transitions between ferromagnetic phases of dif- 
ferent densities remain. The external field deforms the 
phase diagram of a magnetic fluid shifting coexistence 
lines between these phases. Kawasaki studied a magnetic 
lattice gas pj] , which implements one of the ways to model 
properties of magnetic fluids. On temperature-density 
phase diagrams in Ref. jjj one can see that the gas-liquid 
binodal of the magnetic fluid significantly lowers after 
the application of an external magnetic field. Vakarchuk 
with coworkers M and Lado et al. ||^] studied another 
model of magnetic fluids — the fluid of hard spheres with 
embedded Heiscnberg spins. They concluded that in such 
a system the temperature of the gas-liquid critical point 
(the top of binodal) increases after the application of an 
external field. The nature of such a discrepancy can be 
various, because the types of models (continuum fluids 
and the lattice gas 0) as well as approximations 
used differ. 

In this letter results of a more detailed investigation of 
the magnetic (Ising) lattice gas are reported. The first 
new point is the inclusion of the nonmagnetic interaction 
between particles. Another one consists in overcoming 
of limitations applied by the mean field approximation 
(MFA). ft is known that this approximation is good for 
very long-range potentials only (and becomes accurate 
for a family of the infinite-ranged ones, the so-called Kac 
potentials (|||). The MFA can not reproduce some es- 
sential features of the systems with the nearest-neighbor 
interaction (such as the percolation phenomena in the 
quenched diluted Ising model, differences in magnetic 



properties of the quenched site-disordered Ising model 
and the annealed one 0). Such drawbacks can be over- 
come with the two-site cluster approximation (TCA) || . 

We shall show within both the MFA and the TCA 
that at different values of the nonmagnetic interaction 
between particles the Ising fluid demonstrates various ef- 
fects of external fields. In particular, at moderate values 
of the nonmagnetic attraction the field effect on the crit- 
ical temperature is nonmonotoneous. 

In lattice models of a fluid its particles are allowed to 
occupy only those spatial positions which belong to sites 
of a chosen lattice. The configurational integral of a sim- 
ple fluid is in such a way substituted by the partition 
function 

Z = Sp exp (-pH) , [3 = 1/(*bT), 

H = H - fiN = -- } Iijiiirij - n) nj, 



where n iy which equals or 1, is a number of particles 
at site i. Sp means a summation over all occupation 
patterns. The total number N of particles is allowed to 
fluctuate, fi is a chemical potential, which should be de- 
termined from the relation 



TV = 



) H = Z- 1 Sp(---)exp(-f3H)- (1) 



H 



Lattice models, due to particles can not approach closer 
than the lattice spacing allows, automatically preserve 
the essential feature of molecular interaction: nonover- 
lapping of particles. The lattice fluid with nearest neigh- 
bor interaction is known to demonstrate the gas-liquid 
transition only. Nevertheless, the lattice gas with inter- 
acting further neighbors possesses a realistic (that means, 
argon-like) phase diagram with all transitions between 
the gaseous, liquid and solid phases being present 

We shall consider a magnetic fluid in which the parti- 
cles carry Ising spins Si = ±1 and there is also an addi- 
tional exchange interaction between the particles 
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H - fiN = ~ JijSiriiSjnj - h ^ Si: 



i j 

— IijlliUj 



(2) 



Each site can be in three states: empty (i), occupied by 
the particle with spin up (ii) or down (iii). The trace in 
the partition function implies a summation over all 3^ 
states, where Af is a number of sites. 

An interaction of fluctuations are totally neglected in 
the mean field approximation used in the previous studies 
of the model |jj . This can partially be recovered using the 
idea of "clusters" . The partition function of a finite group 
of particles in an external field can be evaluated explicitly. 
A contribution of the other particles may be expressed 
in terms of the effective field, and this field has to be 
evaluated selfconsistently. From such a point of view the 
MFA is a one-site cluster approximation, in which each 
cluster comprises one site. Increasing the size of clusters 
one may expect to obtain more accurate results. Indeed, 
the results of the two-site cluster approximation turns 
to be accurate for the one-dimensional systems jlo| and 
on the Cayley tree 0. Here we shall formulate such an 
approximation for the Ising lattice gas with the nearest 
neighbor interactions. For the sake of brevity we shall not 
use the cluster expansion formulation which has some ad- 
vantages, such as the possibility to calculate corrections 
of a higher order and correlation functions of the model 
||. Instead we shall rely on the first order approxima- 
tion and closely follow the derivation by Vaks and Zein 
0. Let us introduce the effective-field Hamiltonian of a 
single site 



H, 



-hSiUi 



[m t = Ht — —hSiUi - jSnj, (3) 

where h = h + zip, p, = ji-\-zip, (p and ip are effective fields 
substituting for interactions with nearest neighbor sites, 
z is a first coordination number of the lattice. In the 
two-site Hamiltonian the interaction between a pair of 
the nearest neighbor sites is taken into account explicitly 

H% j fJiTli fJiTZj — j — J ij Si Ttx Sj TL j h S%TL'i h S j Tt j 

-Iijmrij - Jxrii - fr'rij, (4) 

where h! = h+z'ip, fi' = fi+z'if), and z' = z — 1 due to one 
of the neighbors is already taken into account. The fields 
have to be found from the selfconsistency conditions that 
require an equality of average values calculated with the 
one-site and two-site Hamiltonians. To determine <p and 
if) it is sufficient to impose these conditions on the aver- 
age values of spin m — (SiUi)^ and of the occupation 
number n = {ni) n , 

(ni) n . = (ni) n .. ; (S l n i ) H[ = (S l n l ) H ^ . (5) 

This approximation leads to the following expression for 
the internal energy 



U/Af = ~-^Jq (SiHiSjUj) 
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h (SiUi) 



Hi 



--IoK«j) H .. , 



(6) 



zV are inte- 



where Jo = J2j Jij = z ^ an d 

gral interaction strengths, K and V denote, respectively, 
the magnetic coupling and the nonmagnetic attraction 
between nearest neighbor sites. The expression (||) can 
be computed explicitly in terms of the fields ip and tp and 
model parameters. The other thermodynamic potentials 
can be found in a straightforward way. For example, the 
grand thermodynamic potential f2 of the model satisfies 
the following Gibbs-Hclmholtz equation 



U-fiN. 



(7) 



The solution of this differential equation, taking into ac- 
count relations (||), reads 

PQ/Af= z'mSpexp(-/3Wi) - | lnSpexp(-/3H«). (8) 

It is possible to build isotherms of the fluid using the 
thermodynamic relation Q = —PV and expression (|^) 
and solving the system of nonlinear selfconsistency equa- 
tions (|^). It is possible and convenient to exclude the 
chemical potential /i and the field ij) from final equations 
of state: 

[3PV/KT = -ln(l -n) + | ln[(l - n)(l - x) + xr], (9) 



tanh/3/i = p 



sinh 20h' 



cosh2/W + exp(-2/3iT) 
+ (1 -p) tanh/3ft/, 



(10) 



where 



p = 1 - (1 - n)/r, r = 0.5 + y/ (n - 0.5) 2 +n(l - nj/x 
2 exp(-/3V - pK) cosh 2 ph' 



cosh2/W + exp(-2/?.?f) 



(11) 



p is a probability that a randomly chosen nearest neigh- 
bor site to a given particle is occupied. In the limit 
z oo and K — > 0, V — > (keeping Jo and Vq con- 
stant) TCA formulae (P-|ll|) turn into the results of the 
MFA. 

At low temperatures the isotherms of the fluid con- 
tain the "liquid" and "gaseous" parts separated by a re- 
gion of the negative compressibility. The thermodynamic 
states, in which the compressibility of the uniform fluid 
is negative (|p < 0), constitute the spinodal region on 
the temperature-density phase diagram. In this region 
the fluid is thermodynamically unstable and must sep- 
arate on phases of different densities. The densities of 
coexisting phases can be determined with the Maxwell 
rule of areas applied to the nonmonotoneous sections of 
the isotherms. Also at sufficiently low temperatures and 
h = an isotherm of the model has a break at the den- 
sity in which nonzero solution of selfconsistency equation 
@ for the field ip appears, and the second order phase 
transition to the ferromagnetic phase occurs. 
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FIG. 1. The phase diagram of the magnetic (Ising) lattice gas within the mean field (a) and two-site cluster (b) approx- 
imations, n = N/N and t = ksT/Jo are reduced density and temperature, respectively. Picture b represents TCA results 
for the model with the nearest neighbor interactions on the simple cubic lattice (z — 6). The results of MFA (picture a) 
are independent of interaction range and lattice structure. Three families of lines are built for three values of nonmagnetic 
interaction v — Io/Jq. The thick lines are binodals and Curie lines at zero external field. The thin lines with attached numbers 
represent the binodals at nonzero values of the external field h/Jo. 



Figure |l|a shows the temperature-density phase dia- 
gram of the model within the mean field approximation 
at different model parameters. There are three fami- 
lies of lines for three values of the nonmagnetic inter- 
action strength v — Iq/Jq. The bold lines correspond 
to zero field case (h = 0). The straight line is a line of 
Curie points that separates paramagnetic and ferromag- 
netic regions. Within the mean field approximation the 
slope of the Curie line is independent of v. Under bin- 
odals (convex lines resting on the points (0,0) and (1,0) 
on the phase diagram) a phase separation takes place: 
one of the phases (vapor) is rarefied, the other (liquid) is 
denser. The nonmagnetic attraction between particles, of 
course, favors the phase separation — the binodal moves 
upward with increasing v. At small v the phase diagram 
possesses of the tricritical point: the top of binodal lays 
on the Curie line, the liquid is ferromagnetic and the gas 
is paramagnetic. The tricritical point disappears at large 
v, in this case the top of binodal deviates from the Curie 
line in the paramagnetic region. 

The external field dissolves the ferromagnetic transi- 
tion and therefore eliminates the Curie line. The gas- 
liquid binodals in presence of the external magnetic field 
are depicted with thin lines. The attached numbers are 
strengths of the field hj Jq. At small v there is a temper- 
ature interval, where the external field suppresses phase 
separation — the top of binodal shifts downward in agree- 
ment with the results of Kawasaki M. Nevertheless, at 
large nonmagnetic attractions (e.g., v = 4) the reverse 
field effect takes place. At moderate values of the non- 
magnetic interaction (for example, v = 2) the field effect 
is nonmonotoneous — weak fields lower the top of bin- 



odal, stronger fields shift it up. 

In Fig. |l|b one can see that the TCA, besides quan- 
titative differences, gives some qualitative corrections to 
the MFA results. Within the TCA the Curie line be- 
comes slightly concave, and the nonmagnetic attraction 
between particles increases the Curie temperature. The 
latter effect can be justified by the qualitative arguments. 
Indeed, the nonmagnetic attraction increases the proba- 
bility that a randomly chosen pair of the nearest neigh- 
bor sites is occupied. Since at this sites particles interact 
magnetically, the magnetic interaction becomes more ef- 
fective, and the Curie temperature increases also. There- 
fore the account of density fluctuations in the TCA leads 
to the dependence of the Curie temperature on v. In 
the case of the non-compressible fluid (n — 1) the den- 
sity fluctuations are absent, and the Curie temperature 
is independent of the nonmagnetic attraction. 

The TCA predictions concerning the effect of field sup- 
port the MFA results. The variety of field effects may 
be explained by the existence of two concurrent ten- 
dencies. The first, the external field aligns the spins, 
which leads to the more effective attraction between par- 
ticles (let us remind that at v = particles with parallel 
spins (SiSj = 1) attract and those with opposite spins 
(SiSj = —1) repulse). This raises the binodal (for ex- 
ample, in simple nonmagnetic fluids the binodal goes up 
when the interaction increases). The second tendency 
takes place, if the susceptibility of the rarefied phase is 
larger than that of the coexisting dense phase. In this 
case the magnetization and, consequently, the effective 
attraction between particles grow better in the rarefied 
phase. This decreases the energetical gain of the phase 



separation. Therefore the second tendency suppresses 
the gas-liquid separation in the fluid and counteracts the 
first tendency. The second tendency is very strong at 
ft = and v = in the region of the tricritical point, 
where the vapor (paramagnetic) branch of binodal almost 
coincides with the Curie line (where the susceptibility 
tends to infinity), whereas the branch of the coexistent 
liquid phase rapidly deviates from the Curie line. As a 
result, the external field lowers the top of the binodal. 
The second tendency gets weak and disappears when the 
susceptibilities of the coexistent phases levels; they are 
comparable, for example, if both liquid and vapor phases 
are paramagnetic. The behavior of v — 4 binodal (see 
Fig. [j]) demonstrates this feature. A relation between the 
susceptibilities results from various factors. For example, 
a short-rangeness of the interactions levels the suscepti- 
bilities and weakens the second tendency Jll[]. It can be 
seen from the following observation of the field effect at 
v = 2: in Fig. [j]a the top of binodal at ft = 0.01 is higher 
than that at ft = 0.1, whereas in Fig. [j]b the reverse sit- 
uation takes place. Since the results of the MFA (as well 
as those of the TCA in the limit z — > oo) are correct 
for the long-ranged potentials, whereas for our case the 
TCA is much more accurate, the corrections provided by 
the TCA have to be attributed to differences between the 
systems with the long-range and short-range potentials. 



Still more illuminating confirmation of the "bi- 
tendency" explanation one can see in Fig. [|. There is 
the phase diagram of special topology, which takes place 
at intermediate v. The model at ft = and t = 0.63 
undergoes two first-order phase transitions. At this tem- 
perature the fluid can be in three phases: paramagnetic 
gas (at n < 0.35), paramagnetic liquid (0.65 < n < 0.69), 
ferromagnetic liquid (n > 0.83). What we would like to 
emphasize is that weak external fields (e.g., ft/ Jo = 0.01) 
raise the binodal at n — 0.5 and lower it at n = 0.75. 
Such behavior completely fit into the "bi-tendency" ex- 
planation: at n = 0.5 and ft = both phases are para- 
magnetic, the second tendency is absent, therefore the 
external field favors the phase separation; at n = 0.75 
the second tendency wins at small fields, like in the case 
v = 2 (see Fig. |Tj). 

One can see that the lattice gas approach may be suc- 
cessfully used for description of complex fluids when con- 
tinual approaches lead to too complex calculations or do 
not give satisfactory results. It is this situation that takes 
place when one determines the influence of the nonmag- 
netic attraction on the Curie temperature [jll],[l2| ■ In this 
case the account of the short-range correlations within 
the cluster approach yields qualitatively new results in 
comparison with current continual methods. 
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FIG. 2. The phase diagram of the model at v = 3 within 
the TCA. The inset shows the vicinity of the top of h — 
binodal in a more detail. 
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